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[1 ] COILS the CFMIP-OASS Intercomparison of Large Eddy Models (LESs) and 
single column models (SCMs)investi gates the mechanisms of cloud feedback in 
SCMs and LESs under idealized climate change perturbation. This paper describes the 
COILS results from 15 SCMs and 8 LES models. Three cloud regimes over the sub- 
tropical oceans are studied: shallow cumulus cumulus under stratocwnulus, and well- 
mixed coastal stratus/stratocumulus. In the stratocumulus and coastal stratus regimes, 

SCMs without activated shallow convection generally simulated negative cloud feed- 
backs while models with active shallow convection generally simulated positive cloud 
feedbacks. In the shallow cumulus alone regime this relationship is less clear likely 
d ue to the changes in cloud depth lateral mixing and precipitation or a combination 
of them. The majority of LES models simulated negative cloud feedback in the well- 
mixed coastal stratus/stratocumulus regime, and positive feedback in the shallow 
cumulus and stratocwnulus regime. A general framework is provided to interpret 
SCM results: in a warmer climate, the moistening rate of the cloudy layer associated 
with the surface-based turbulence parameterization is enhanced; together with weaker 
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large-scale subsidence, it causes negative cloud feedback. In contrast, in the warmer cli- 
mate, the drying rate associated with the shallow convection scheme is enhanced . This 
causes positive cloud feedback. These mechanisms are summarized as the " ESTS" 
negative cloud feedback and the SCOPE” positive cloud feedback (Negative feedback 
from Surface Turbulence under weaker Subsidence-Shallow Convection PositivE 
feedback) with the net cloud feedback depending on how the two opposing effects 
counteract each other. The LES results are consistent with these interpretations. 

Citation: Zhang, M., et al. (2013), CGILS: Results from the first phase of an international project to mlderstand the physical 
mechanjsms of low cloud feedbacks in single column models, J. Adv. Model. Earth Syst., 5, doi: 10.1 002/20 13MS000246. 


1. Introduction 

[2] Cloud-climate feed backs in General Ci rcula tion 
Models (GCMs) have been the subject of intensive 
stud y for the last four decades [e.g., Randall et al, 
2007]. These feedbacks were iden tified to be one of the 
most significant uncertain ties in projecting future global 
warming in past IPCC (Inter- Governmental Panel for 
Climate Change) Assessment Reports (AR), as well as 
in coupled model sim ulations that will be used for the 
upcoming AR5 [ Andrews et al, 2012]. Despite much 
progress toward understanding cloud feed backs [ Bony 
et al, 2006], however, there is still a general lack of 
knowl edge a bou t thei r mechanisms . Understanding the 
physical mechanisms is necessary to increase our confi- 
dence in the sensitivity estimates of climate models. 

[3] Cloud-dinrnte feed backs refer to the radiative 
inlpact of changes of clouds on climate change. Because 
clouds are not explicitly resolved in GCMs, they are the 
product of an interactive and elaborate sui te of physical 
parameterization s. As a result, i t has been a challenge 
to decipher cloud feed back mechani sm s in climate mod- 
els. Clouds also interact wi th the resolved-scale a tmos- 
pheric dynamical circulations through their impact on 
latent and radiative heating. 

[4] I n view of the challenges, CLMIP (the Cloud F eed- 
back Model Intercomparison Project) and GASS (Global 
Atmospheric System Studies) initiated a joint project- 
CGILS (the CFMIP-GASS Intercomparison of Large 
Eddy Models (LESs) and single column models (SCMs)) 
to analyze the physical mechanisms of cloud feedbacks in 
SCMs by using an idealized experimental setup. The focus 
of CGILS is on low clouds in the subtropics, because sev- 
eral studies have demonstrated that these clouds contrib- 
ute significantly to cloud feedback differences in models 
[G.g.,Bony and Dufr esne, 2005; Zelinka et al, 2012]. The 
role played by these clouds is consistent with the fact that 
low clouds have the largest net cloud -radiative effect, in 
contrast to deep clouds in which the positive longwave 
and negative shortwave cloud effects largely cancel out 
[G.g.,Rammlathan etal., 1989]. 

[s] The objective of this paper is to describe the 
CGILS project and results from 15 SCMs and 8 LES 
models. Section 2 briefly descri bes the experimental 
design and large-scale forcing data. Section 3 gives a 
brief description of the partici pating model s. Section 4 
discusses simulated clouds and the associated physical 
processes. Section 5 presents cloud feedback results . A 
brief sunlmary is given in Section 6. 


2. Experimental Design and 
Large-Scale Forcing Data 

2.1. Experimental Design 

[6] The CGILS experimental design was described in 
Zhang et al [2012], which is schematically shown in Fig- 
ure 1 . In t he control climate (CTL), sea surface tern per- 
at ure (SST) is specified along the GCSS/WGNE Pacific 
Cross Section Intercomparison (GPCI) [Teixeira et al, 
2011] in the northeast Pacific by using the ECMWF 
(European Center for Medium-Range Weather Fore- 
casts) Interim Reanalysis (ER A-Interim) [ Dee et al, 
201 1] July 2003 condition as given in Table I of Zhang 
et al [2012]. In the perturbed climate, SST is uniformJ y 
raised everywhere by 2° as in Cess et al [1990]. Large- 
scale horizontal advection and vertical motion , corre- 
sponding to the underlying SST, were derived and used 
to force SCMs and LES models. The perturbed climate 
is referred to as P2S, wi th "S" denotes that t he large- 
scale subsidence is also differen t from CTL [Bretherton 
et al, 2013]. The models simulate changes of clouds in 
response to changes of SST and the associated large- 
scale atmospheric conditions. 



Pertnrbed 

SST ^ m 

Figure 1. Schematics of the experimental setup . The 
atmospheric tern perature and wa ter vapor are con- 
structed based on moist adiabat and fi xed relative 
h umidity, respecti vel y . The large-scale subsidence is cal- 
culated based on the clear-sky themlodynamic eq ua- 
tion. These fields change wi th SST warming of 2°C in 
the perturbed climate. 
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Table 1 . Participating Models, Main References, and Contributorsa 






Layers: Total/ 

Models Acronyms 

Model Institution 

References 

Contributors 

(p>700hPa) 

SCM ( 15) 

ACCESS (Australian 

Australian Commonweal th 

Hewittel al. [201 1J 

Charmaine Franklin 

38/12 

Community Climate and 

Scientific and Industrial 




Earth System Simulator) 

Research Organisation/ 
Centre for A ustralian 
Weather and Climate 
Research 




CA M4 (Communiry 

National Center for 

Neale el al [20IOJ 

Minghua Zhang, Cecile 

26/5 

Atmospheric Model 

Atmospheric Research 


Hannay, and Philip 


Version 4) 

(NCAR), USA 


Rasch 


CAMS (Community 

National Center for 

Neale el al. [2012J 

Minghua Zhang, Cecile 

30/9 

Atmospheric Model 

Atmospheric Research 


Hannay, and Philip 


Version 4) 

(NCAR), USA 


Rasch 


CCC (Canadian Centre for 

Canadian Centre for 

Ma et al. [2010] 

Phillip Austi n and Knut von 

35/14 

Climate) 

Climate Modell ing and 
Analysis, Canada 


Salzen 


CLUBB (Cloud Layers 

University of Wisconsin at 

Goto= et of. [2002a, 2002J, 

Vincent Larson and Ryan 

41129 

Unified By Binormals) 

Milwaukee, USA 

LJHson and Galaz [2005], 
and Galaz el al. [2007] 

Senkbeil 


ECHAM6 (ECMWF- 

Max-Planclc Institute of 

Ste liens el al. [2013 J 

Suvarchal Cheedela and 

31/9 

University of Ham burg 
Model Version 6) 

Meteorology, Germany 


Bjorn Stevens 


ECMWF (European Center 

European Center for 

Neggers et al. [2009a, 2009bJ 

Martin Koehler 

91120 

for Medium Range 

Medium Range Weather 




Weather Forecasting) 

Forecasti ng 




EC-ETH (ECMWF- 

Swiss Federal I nstitute of 

ll-otta et al. [201 1J 

Colombe Siegenthaler-Le 

31/9 

Eidgenossische 

Technology , Switzerland 


Drian, Isotta Francesco 


Technische Hochschule) 



Alessandro, and Ulrike 
Lohman 


GFDL-AM3 (Geophysical 

NOAA Geophysical Fluid 

Donner et al. [20 1 1 J 

Jean-Christophe Golaz and 

48/12 

Fluid Dynamics 

Dynamics Laboratory, 


Ming Zhao 


Laboratory Atmospheric 
Model Version 3) 

USA 




GISS (Goddard Institute for 

NASA Goddard Institute 

Schmidt et al. [2006J 

Anthony DelGenio and 

40/9 

Space Snldies) 

for Space Sn.dies, USA 


Audrey Wol f 


GMAO (NASA Global 

NASA Goddard Space 

Rienecker et al. [2008J and 

Andrea Molod , Max 

72/13 

Modeling and 

Flight Center, USA 

Mofod et af. (2012] 

Suarez, and Julio 


Assimilation Office) 



Bacmeister 


HadGEM2 (Hadley Centre 

Met Office, United 

Lock et al. [2001] and 

Adrian Lock and Mark 

38112 

Global Environment 
Model version 2) 

Kin gdom 

Mar tin et al. [20 1 1 J 

Webb 


JMA (Japan Meteorological 

Japan Meteorological 

Kawai [2012J 

Hideaki Kawai 

60/16 

Agency) 

Agency , Japan 




1PSL (Institute Pierre Simon 

Institute Pierre Simon 

Hourdin et al. [2006J 

Florent Brient, Sandrine 

39/12 

Laplace) 

Laplace (IPSL), France 


Bony, and Jean-Louis 
Dufresne 


RA CMO ( Regional 

Royal Netherlands 

Neggers et al. [2009a, 2009bJ 

Roel Neggers and Pier 

91120 

Atmospheric Climate 

Meteorological Institute, 


Siebesma 


Model) 

the Netherlands 




LES (8) 
DALES (Dutch 

R oyal Netherlands 

Heus et al. [20 1 OJ 

Stephan de Roode 


A tmospheric Large-Edd y 

Meteorological Institute, 




Simulation) 

the Netherlands 




LARC (NASA Langley 

NASA Langley Research 

Xu et al. [20IOJ 

Anning Cheng and Kuan- 


Research Center) 

Center, USA 


man Xu 


SAM (System for 

University of Washington/ 

Klwirtol ltdino 11 and Randaff 

Peter Blossey, Chris 


Atmospheric Models) 

Stony Brook University, 

[2003J 

Bretherton, and Marat 



USA 


Khairoutdino v 


SAMA (System for 

University of Washington/ 

Khairtol ltdino 11 and Randaff 

Peter Blossey, Chris 


Atmospheric Models) 

Stony Brook University, 

[2003J and Bfossey et al. 

Bretherton, and Marat 



USA 

[201 3J 

Khairoutdinov 


MOL EM (Met Office Large 

Met Office, United 

Lock [2009J 

Adrian Lock 


Eddy Model) 

Kingdom 




MOLEMA (Met Office 

Met Office, United 

Lock [2009J and Bfoss ey 

Adrian Lock 


Large Eddy Model) 

Kingdom 

etal. [2013] 



UCLA (University of 

Max Plank Institute of 

Stellens et al. [2005J and 

Thijs Heus, Irina Sandu, 


California at Los 

Meteorology , Germany/ 

Stevens and Seifert [2008 J 

and Bjorn Stevens 


Angeles) 

University of California 
at Los Angeles, USA 
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Table 1. (continued) 

Model s Acronyms 

Model Institution 

References 

Contributor s 

Layers: Total/ 
(p > 700 hPa) 

WRF (Weather Research 
and Forecasting) 

National Center for 
Atmospheri c Research/ 
Brookhaven National 
Laboratory 

Endo et al. [201 1) 

Satosh End and Yangang 
Liu 



"The number of vertical layers and layers between the surface and 700 hPa for SCMs are given in the last column. 


[7] Three locations along the GPCI cross section are 
selected for study. They are labeled as S6, SI 1, and S12 in 
Figure 2, which also shows the distribution of low cloud 
amount in the sununer (JJA, June to August) from t he 
merged CALIPSO, CJ oudSat, CER ES, and MODIS sat- 
ellite prod uct C3M [ Kato et al., 2011 ; Xu and Cheng, 
2013]. Typical regimes of clouds at these three locations 
are shallow cumulus (S6), cumulus under stratocumulus 
(SI 1), and well-mixed stratocumulus or coastal stratus 
(SI 2). On the basis of dominant cloud types, they are 
referred to as shallow cum ul us, stratocumul us, and 
coastal strat us, respectively. The locations and val ues of 
summer-time surface meteorological variables in the con- 
trol climate can be found in Table 1 of Zhang et al. [2012]. 

2 . 2 . Forcing Data 

[s] The SCM and LES forcing data refer to the large- 
scale horizontal advective tendencies and vertical veloc- 
ity, and surface bound ary condi tions that are specified 
in the model sim ulations. The SCMs calculate the time 
evol u tion of water vapor and temperature as follows 


l Randall and Cripe , 1999]: 



• ■'(■ 



5 2 Ji” A7h 

2x (gJ, 11 - 

(i) 

@ @ 

w w 1/iy :(. 

LS LS @ ' 

(2)111 5 @111 2 Ji’ . \7q 

2 x ® m - 

(2) 

@ @ pHy 

LS LS @ • 


where h and q are potential temperat ure and water vapor 
mixing ratio . Subscript '^"denotes model calculations; 
"LS" stands for large-scale; other symbols are as com- 
monly used. The Erst teml on the right-hand side (RHS) 
of equations ( I) and (2) is calculated from physical 
parameterizations (with subscript "phy s"). The last two 
terms contain the specified large-scale horizontal advec- 
tive forcing and subsidence. In LES models, conservative 
variables like liquid water potential temperature and 
total liquid water are typically used as prognostic fields 
[e.g., Siebesma et al., 2004; Stevens et al., 2005]. Equa- 
tions (1) and (2) represent domain averages. The atmos- 
pheric winds and initial relative humidity are specified 
by using the ERA-Interim for July 2003. Initial profiles 
of atmospheric temperature are assumed to follow moist 
adia bat over the warm pool and weak gradient approxi- 
mations at other locations [Sobel et al., 200 I]. Surface 
latent and sensi ble heat fluxes are calculated internally 
by each model from the specified SST and winds. 

[9] The large-scale horizontal advective tendencies 
and subsidence in equations (1) and (2) are specified 


according to SST. In the free troposphere , they are 
deri ved based on the clear-sky thermodynamic and 
water vapor mass continuity equations , in which radia- 
tive cooling in the thermod ynamic equation is balanced 
by subsidence warming and horizontal ad vection, with 
the radiati ve cooling calculated by usi ng the RRTM 
radiation code [ Mlm ver et al., 1997] and the horizontal 
advection constrained by ERA-Interim. Below the alti- 
tude of 900 hPa, the horizontal advecti ve forcing of 
tern perature and water vapor are calcu lated using the 
SST spatial gradient and specified surface relative 
h umidity. The detailed derivation of the CGILS forcing 
data and its comparison with the corresponding GCM 
and ERA-Interim can be found in Zhang et al. [2012]. 

[10] Figure 3 a shows the derived vertical profiles of 
x LS in CGILS CTL (solid lines) and ERA -Interim 
(dashed lines) at the three chosen locations. The 
obtained values match well with ER A-Interim in the 
lower troposphere. Among the three locations, the subsi- 
dence rate is the strongest at S12 and the weakest at S6. 

[11] Figure 3b shows the comparison of the derived 
x LS between CTL (solid lines) and P2S (dashed lines) 
used in the sim ulations. It is seen that subsidence is 
weaker in the warmer climate. Figures 3c and 3d show 
the corresponding profiles of horizontal advective ten- 
dencies of temperature and water vapor, respectively. 
In the free troposphere, these profiles , along with 



°170W 160 W 150W 140W 130W 120W HOW 


10 15 20 25 30 35 40 45 50 SS 60 65 70 75 80 

Figure 2. Averaged amount of low clouds in June- 
July-August (%) from the C3M satellite data. The red 
line is the northern por tion of the GPCI (see text); the 
symbols ”S6," "Sll ," and "S12" are the three locations 
studied in the paper. 
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Figure 3. (a) Large-scale pressure vertical velocity at the three locations in the con trol dinmte (solid lines), and in 
the ER A-Interim (dashed), (b) Same as Figure 3a except that the dashed 1 ines denote subsidence rates in the 
warmer climate, (c) Same as Figure 3b except for horizontal advective tendency of temperature, (d) Same as Figure 
3c except for advective tendency of water vapor. 


the profiles of X u;, SST, and initial atmospheric tem- 
perature and water vapor, satisfy the clear-sky atmos- 
pheric thermod ynamic and water vapor mass continuity 
equations under 15 July insolation conditions. Zhang 
et al. [2012] showed that the changes in the forci ng data 
between CTL and P2S in Figure 3 capture the essential 
features in GCMs . All data are available at the 
CGILS websi te http ://atmgcm.msrc.sunysb .edu/cfmip_ 
figs/ Casespeci fication .htm I . 

2.3. Simulations 

[12] We use the change of cloud-radiative effect 
(CRE) from CTL to P2S, as in many previous studies, 


to measure cloud feed backs. Even though Soden et al 
[2004] suggested other better diagnostics of cloud feed- 
backs, CR E is used for simplicity, which should not 
affect the results of this paper. 

[13] The SCMs and LES are integrated to q uasi- 
equilibrium states by using the same steady large-scale 
advective tendencies and subsidence as forcing data. 
Each model ran six simulations: CTL and P2S at the 
three locations of S6, SI 1 , and S 12. Since the forcing is 
fixed, a model may event ually d rift if its radiative cool- 
ing rate in the free a tmosphere d iffers from the rate 
used in the derivation of the prescribed large-scale sub- 
sidence. To prevent models from similar drifting, at 
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Table 2. Boundary-Layer Tw-bulence Schemes in SCMs 


Models 

R eferences 

Local 

Kc 

Cloud-top 
Entrai nment 

Cotmter 

Gradient 

C,. 

ACCESS 

Lock etal. [2000] 

N 

y 

y 

CAM4 

Holtslag and Bal *ille [1993] 

N 

N 

y 

CAMS 

Br etherton and Park [2009] 

y 

y 

N 

CCC 

P0J1 Sal en et al. [20 13] 

y 

y 

y 

CLUBB 

Galazet al. [2002a, 2002], Larson and Galaz [2005], 
and Galaz et al. [2007] 

N 

N 

N 

ECHAM6 

Ste Pens etal. [2012) 

y 

N 

N 

ECMWF 

Neggers et al. [2009a, 2009b] and Lock [2000] 

N 

y 

y 

EC-ETH 

Brinkop and Roeckner [1995] 

y 

N 

N 

GFDL-AM3 

Lock etal. [2000) and Louis and Geleyn [1982) 

N 

y 

N 

GISS 

Holtslag and Moeng [1991) and Del Genia et al. [1996] 

y 

y 

y 

GMAO 

Lock etal. [2000) and Louil- and Geleyn [1982) 

N 

y 

y 

HadGEM2 

Lock etal. [2000] 

N 

y 

y 

JMA 

Mel/or and Yamada [ 1974] and Kalllai [2012] 

y 

N 

N 

rPSL 

Hourdin et al. [2006] 

y 

N 

y 

RACMO 

Neggers et al. [2009a, 2009b] 

N 

y 

y 


pressure less t han 600 hPa , temperatu re and water 
vapor mixing ratio are relaxed to thei r initial conditions 
wi th a time scale of 3 h. In LES models , they are relaxed 
a t altitudes above 4000 m for S6, 2500 m for SI 1, and 
1200 m for SI 2, respectively, to reduce computational 
costs and allow for high vertical resolu tions in shallow 
domains. Some LES models did not com plete all six 
sim ulations. 

[1 4] Most of the SCMs are integrated for 100 days. 
Based on a visual inspection of statistical eq uili brium, 
the averages of their last period of about 50 days are 
used . Most LES simulations reached quasi-equilibrium 
sta tes after 10 days, in which case the last 2 d ays are 
used in the analysis. Zhang and Br etherton [2008] ana- 
lyzed the transient behavior of the Comm uni ty Atmos- 
pheric Model (CAM) under constant forcing and 
showed that the interaction of djfferent physical param- 
eterization com ponents can create quasi-periodic 
behavior s of model sinm lation wi th time scales longer 
than a day. Since LES models contain fewer parameter- 
ization com ponents, the ini pact of this type of interac- 
tions is reduced , wruch may explain why LES models 
reach quasi steady states in shorter time than SCMs. To 
our knowled ge, CGILS is the first LES intercom parion 
study to investigate clouds by integrating them to 
quasi-equilibrium states. 

3. Models and Differences in Physical 
Parameterizatio ns 

[1 s] Fifteen SCMs and eight LES models pa r ticipated 
in this study. Many parent GCMs of the SCMs also 
par ticipated in the Coupled Model Intercomparison 
Project 5 (CMIP5) . Table 1 lists the model names , main 
references, and CGILS contri bu tors. It also gives the 
num ber of total vertical model layer s and n um ber of 
layers between the surface and 700 hPa in SCMs. The 
SCM vertical resol ution in the boundary layer (PBL) is 
generally not sufficient to resolve the observed or LES 
simulated thin stratocumulus clouds. No attempt is 


made to make them finer since our objective is to under- 
stand the behavior of operational GCMs. For the LES 
mod els, however, because they are intended as bench- 
marks, much higher resolutions are used . The horizon- 
tal resolutions of LES models are 100 m, 50 m, and 25 
m, respecti vely, at S6, S 11, and S 12. The vertical resolu- 
tions of the majori ty of LES are 40 m, 5 m, and 5 m, 
respectively, at the th ree locations. More detailed 
descriptions of the CGILS LES models are given in a 
companion paper by Blossey et al. [2013]. 

[16] The physical parameterizations in the SCMs rele- 
vant to the present stud y are the PBL, shallow convec- 
tion, and cloud schemes. For PBL schemes, the generic 
form can be wri tten in terms of t urbulent flux a t the 
mod el interfaces : 

♦ . ♦ 

t»’ S'5 2 Ke g2c ( . ; (3) 

where z is height , w is vertical velocity, S is a conserva- 
tive model prognostic variable. Prime represents the 
turbulent perturbation from the mean that is denoted 
by the overbar. Kc is the eddy diffusivity, and CC is the 
counter-gradient transport teml. In addition to resolu- 
tion, the differences in PBL schemes among the models 
are in their formulations of Kc and CC. For Kc, some 
mod els parameterize i t by using local variables a t the 
resolved scales, such as local Richardson num ber in the 
so-called fi rst order closure models, or local turbulent 
eddy kinetic energy (TKE) [Mellor and Yamada, 1974]. 
Other models use nonlocal em pirical parameterization 
of K c as a function of heigh t relative to the boundary 
layer depth. Another Kc difference among the models is 
its parameterization at the top of the PBL. While some 
mod els have explicit parameterizations of turbulent 
entrainment based on parameters such as cloud -top 
radiative and evaporative cooling, others do not con- 
sider entrainment. For the counter-gradient teml cc, 
some models calculate i t based on surface buoyancy 
fluxes, while others do not have this term. Table 2 
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Table 3. Shallow Convection Schemesa 


Models 

Acronyms 

References 

Trigger 

Lateral 

Entrainment 

Lateral 
Detrainmen t 

Closure 

ACCESS 

Gregory and Rollllltree 
[ 1990) and Grant [200 1) 

Undiluted parcel 

Specified 

Specified 

TKE 

CAM4 

Hack [1994] 

Undiluted parcel 

N 

N 

CAPE 

CAMS 

Park and Bretherton [2009] 

CIN 1TKE 

Buoyancy sorting 

Buoyancy sorting 

CIN 1TKE 

CCC 

iio/j Salzen et al. [2012], 
van Salzen and McFarlane 
[2002], and Grant [2001) 

Undiluted parcel 

Buoyancy profile 

Buoyancy profile 

TKE 

CLUBB 

Galaz et al. [2002a, 2002), 
Larson and Galaz [2005), 
and Galaz el al. [2007) 

N 

N 

N 

High-order 

bi-normal 

distribution 

ECHAM6 

Tiedtke [1989) 

Diluted parcel 

Specified 

Specified 

Moisture con 
vergence 

ECMWF 

Tiedtke [1989) 

Diluted parcel 

Specified 

Diagnosed 

Subcloud moist 
static energy 

EC-ETH 

Von Salzen and McFarlane 
[2002), 

Grant [2001), and Isotta 
el al. [201 1) 

Undiluted 

Buoyancy profile 

Buoyancy profile 

TKE 

GFDL-AM3 

Breth erton and Park [2009) 
and Zhao et al. [2009) 

CIN 1 TKE 

Buoyancy sorting 

Buoyancy sorting 

CTN 1TKE 

GISS 

Del Genia and Yao [ 1993) and 
Del Genia et al. [2007) 

Undiluted parcel 

Buoyancy and speed 

Above neutral level 

Cloud-base 

buoyancy 

GMAO 

Moorthi and Suarez [1992] 

Undiluted 

Diagnosed 

N 

CAPE 

HadGEM2 

Gregory and Roivnlree [1990) 
and Grant [200 1] 

Undiluted parcel 

Specified 

Specified 

TKE 

JMA 

Pan and Randall [1998) 

Diluted parcel 

Di agnosed 

N 

Prognostic 

IPSL 

Emanuel[ 1991 , 1993] 

Undiluted parcel 

Buoyancy sorting 

Buoyancy sorting 

CAPE 

RACMO 

Neggers et al. [2009a, 2009b] 

Unified with PBL 
scheme 

Unified with PBL 
scheme 

Unified with PBL 
scheme 

Unified with 
PBL scheme 


"Some models use the same schemes for deep convections. 


categorizes the PBL schemes in the SCMs accord ing to 
the above at tributes. Cloud-top entrainment in Table 2 
refers to explicit parameterization . PBL schemes form u- 
lated using moist conserved variable and TKE closure 
(such as ECHAM6) may implicitly contain cloud-top 
entrainment. As can be seen, a wide variety of PBL 
parameterizations are used in the SCMs. Because of 
coarse vertical resolutions , however , some of these dif- 
ferences d o not make as m uch of an impact on cloud 
sim ulations as they would if higher vertical resol utions 
were used. 

[11] The majority of SCMs used mass-flux shallow 
convection schemes. The generic foml of convective 
transport for a conservative variable q 1 in these schemes is 


••V 5 M(z)(q lc 2 q, e ); (4) 

where the prime denotes deviation of the bul k proper- 
ties of clouds from the mean ; M is the convective mass 
flux; subscripts c and e represent values in the parame- 
terized cloud model and in the environment air, respec- 
tively. The convective mass flux is calculated from 
parameterized rates of entrainment and detrainment d 

@.7 5 k2 cl: 

M @ 

[Is] Some models do not separatel y parameterize 
shallow and deep convection. The schemes can differ in 
their entrainment and detrainment rates, the closure 


that determines the amount of cloud base mass flu x, 
and convection triggering condition as well as origina- 
tion level of convection. Table 3 categorizes the convec- 
tive schemes in the SCMs based on these mai n 
attrib utes. Among the SCMs, CLUBB, and R ACMO 
use a single scheme to parameterize PBL turbulence 
and shallow convection. 

[19] Cloud schemes in SCMs include a macrophysical 
and a microphysical com ponent. Cloud macro physica 1 
schemes parameterize cloud amount and the grid-scale 
rate of condensation and evaporation. These schemes 
can be generically described by assuming that the total 
water in the air, q v obeys a probability distribution 
function (pdl) P ( qi) within a model grid box. The cloud 
amount is then 

CS P(ql)dql ; (5) 

where qs is the saturation vapor pressure at cloud tem- 
perature. Cloud liquid water qi is then 

-DO 

q,5 (ql2 qs) P(ql )dql : (6) 

[20] Therefore, cloud fraction and cloud liquid water 
are often propor tional to each other in individ ual mod- 
els when the cloud fraction is less than 100%. The cloud 
microphysics scheme treats how condensed water is 
converted to precipitation. Inmost parameterizations, 
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Figure 4. (a-c) Averaged profiles of cloud amount (%) by SCMs for S6, SI 1 , and SI 2, respectively (from top to 
bottom plots), (d-f) Same as Figures 4a-4c but by the LES models, (g-i) From the C3M satellite measurements. 
The blue lines are ensemble averages; the red lines are the 25% and 75% percentiles. 


preci pitation is typically proportional to cloud water, 
which is fur ther proportional to rate of large-scale 
condensation. 

4. Simulated Clouds and Associated Physical 
Processes 

[21] Before investigati ng cloud feed backs, we first 
examine the sim ulated clouds in CTL. Figure 4 shows 
the time-averaged cloud profiles in all 15 SCMs and all 
LES models, with the shallow cumulus location S6 in 
the top row and the stratus location S 1 2 in the bottom 
row. SCMs resul ts are in the left column; LES models 
in the middle col umn; observa tions from C3M for the 
summers of 2006-2009 in the righ t col umn. Note that 
the observations may have categorized drizzle as clouds, 
therefore having a different definition of clouds from 


that i n the models. The bl ue lines denote the ensem ble 
averages or m ultiyear averages; the red lines denote the 
25 and 75 percentiles. Figu re 5 shows exam pies of the 
time -pressure cross sections of these cloud amount from 
a sample of three SCMs (JAM, CAM4, and GISS), 
which are selected because they span the range of model 
differences as will be shown later, and from one LES 
(SAMA). 

[22] Despite large differences among the models, the 
relative rank of cloud-top height and cloud amount a t 
the three locations is correct. The spread in the LES 
mod els is much smaller than that among the SCMs. At 
S 1 1 , LES models sim ula ted cumulus under stratocum u- 
1 us. The use of the steady forcing for all models may 
have amplified the intermodel differences, since in both 
GCMs and the real atmosphere the large-scale circula- 
tion can respond to local differences in the inversion 
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Figure 5. Examples of time evolution of cloud amount (%) simulated by JMA (left column) for S6, SI 1, and S12, 
respectively, from top to bottom plots; CAM4 (middle column); GISS (third column); SAMA (right column). 


height by partially compensating them [Blossey et al., 
2009; Bretherton et al., 2013]. 

[23] We find it instructi ve to use the following mois- 
ture budget equation to probe the physical parameter- 
izations responsible for the simulated clouds in the 
SCMs. It is written as: 

• • • • 

@tv 5 @ v I @_ v 2(2) 

@ @ 

'.I. @ C e strtl 

Irtrb t couv I o (7) 

2 ?' • "ilq l X LS @v ; 

LS 0 

where the varia bles are as commonly used, and the 
tend ency terms have been separated into three physical 
terms representing parameterizations of PBL turbu- 
lence ( turb), convection ( conv), large-scale stratiform 
net condensation (c-e), plus the three -d imensional large- 
scale forcing. As wi 11 be shown later, the separa- tion of 
the physical tendency terms hel ps to provide a 
framework of interpreting cloud feed back behaviors in 
the models. We show the three selected models in Fig- 
ure 6 of the time-averaged profiles of these three terms 
at SI 1 in CTL by using the colored solid lines. The 
black 1 ines are t he simulated grid -box mean clo ud liq uid 
water. The solid dots on top of the black lines donate 
the midpoint of model layer. 

[2A] In the JMA model, only two physical terms are 
active (Figure 6a) in addition to the large-scale dynamic 


forcing . The PBL scheme moistens t he bound ary layer; 
the large-scale condensation d ries it. The resid ual is bal- 
anced by the drying from the large-scale forcing. The 
peak altitudes of the "turb" and "c-e" are the same as 
that of the cloud liquid water. Since the PBL scheme is 
always active, the stratiform condensation scheme 
responds to the PBL scheme. In CAM4, Figure 6b 
shows that shallow convection is active in addition to 
the '"turb" and the '"c-e" terms. The shallow convective 
scheme transports the moisture from the boundary 

layer to the free troposphere. In the GISS model , Figure 
6c shows that shallow convection is also active, bu t 
unlike CAM4, the maximum drying of the '" conv" 
is at the same level as the maxim um level of m turb, " in 
the middle of the cloud layer. These differences will be 
shown later as causes of different cloud feed backs in the 
mod els. In Figure 6, the stratiform condensation tern 1 is 
the direct source of cloud water. 

[25] The in termodel differences in Figure 6 are exam- 
ples of how different parameterization assumptions can 
affect the balance of the physical processes and associ- 
ated clouds. The JMA model used the prognostic 
Arakawa- Schubert convection scheme [Pan and Ran- 
dall, 1998] with fixed cloud base level near 900 hPa in 
the model [.!MA , 2013]. As a resul t , convection is not 
active in this case. CAM4 and GISS both used positive 
Convective Available Potential Energy (CAPE) of undi- 
luted air parcels as cri teria of convection. As a resul t, 
shallow convection is more easily triggered in these two 
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Eigure 6. Solid lines are physical tendencies of water vapor (g/kg/day) in three SCMs at SI 1 for the control cli- 
mate, "turf for turbulence scheme, "conv" for convection scheme, "( c-e)" for net large-scale condensation, and 
"qf for the grid-box cloud liquid water (0.1 g/kg) . The black dots show the midpoint of model layers . The dashed 
lines show the corresponding values in the warmer climate, (a) JMA, (b) CAM4, and (c) GISS. 


models . Nevertheless, the assumptions in their shallow 
convection parameterizations are different. For exam- 
ple, CA M4 does not incl ud e lateral entrainment into 
the convective plumes [Hack, 1994], while GISS has lat- 
eral entrainment [Del Genio and Yao, 1993]. 

[26] Tables 4-6 show the simulated surface sensible 
and latent fluxes, precipitation , cloud water path , and 
cloud-radiative effects in the SCMs at SI 2, Sll , and S6, 
respecti vely, in the control climate. Total cloud amount 
is not incl uded in the ta ble since in some models it is 
contaminated by unreali Stic optically thin clouds in the 
upper troposphere. The expected increase of surface 
latent heat fluxes from S 12 to SI 1 and S6 is simulated 
in most models. However , consistent wi th wha t has 
been shown in the vertical profiles of clouds in Figure 4, 


the models differ greatly in their cloud liquid water 
path, and as a result, in the shortwave cloud radiation 
effect. At S 12, some models did not simulate clouds. As 
shown in Zhang et al. (2012] for the GFDL model, this 
unrealistic behavior is related to the use of steady forc- 
ing. When compared with the LES results of Tables 3-5 
in Blossey et al. (2013], the SCM surface latent heat 
fluxes are generally smaller than in the LES models. 
This is likely related to the use of the steady forcing or 
insufficient entrainment mixing in the SCMs. The pre- 
cipitations and the cloud liquid paths in the SCMs span 
a wide range that brackets the corresponding range in 
the LES models. Since the objective of CGILS is to 
investigate the cloud feedback or the response of the 
cloud fields to a warmer climate, we only use Figure 6 


Table 4. Simulated ields in Control Climate and Their Changes in the Perturbed at S12 in SCMsa 


ModellD 

SH 

LH 

PREC 

TGLWP 

SWCRF 

CRE 

ACCESS 

13.8(2 5.8) 

58.9 (2 2.8) 

0.00 (0.00) 

14.2 (2 5.4) 

2 79.4(35.4) 

2 72.2 (32.3) 

CAM4 

24.7(2 0.6) 

48.3 (4.6) 

0.00 (0.00) 

199.4 ( 1 1.0) 

2 210.4 (2 0.6) 

2 215.5 (2 1.0) 

CAMS 

2 6.0 (0.2) 

2.9 (0.3) 

0.00 (0.00) 

0.0 (0.0) 

0.0 (0.0) 

0.0 (0.0) 

CCC 

26.6 (2 3.6) 

54.4(13.1) 

0.51 (2 0.14) 

186.2 (2 82.5) 

2 100.4(17.2) 

2 100.3(19.5) 

CLUBB 

25.8 (2 1 .6) 

64.7 (11.4) 

0.00 (2 0.00) 

77.8(24.2) 

2 176.2(2 18.2) 

2 170.5 (2 18.0) 

ECHAM6 

2 22.8 (1.9) 

62.2 (2.9) 

1.10(0.10) 

98.1 (0.9) 

2 121.4(0.8) 

2 124.1 (1.6) 

ECMWF 

10.1 (2 3.7) 

68.1 (15.4) 

0.00 (2 0.00) 

12.5 (3.8) 

9.9 (2 5.4) 

12.8 (2 4.2) 

ECETH* 

2 27.9 (43.7) 

1.5 (32.8) 

0.00 (0.0) 

0.0 (0.0) 

0.0 (0.0) 

0.0 (0.0) 

GFDLAM3 

2 4.8 (1.1) 

18.9(2.6) 

0.00 (0.00) 

0.0 (0.0) 

0.0 (0.0) 

0.0 (0.0) 

GISS 

11.3 (2 0.5) 

59.9(10.7) 

0.35 (0.22) 

140.9(95.1) 

2 109.0(2 24.5) 

2 108.2 (2 24.3) 

GMAO 

1 .3 (0.2) 

35.5 (2.1) 

0.50 (2 0.50) 

0.8 (2 0.8) 

2 LO (0.9) 

2 1.3 (1.2) 

HadGEM2 

17.0(2 1.8) 

61.2 (7.2) 

0.70 (2 0.30) 

23.9 (2 4.4) 

2 95.5 (13.5) 

2 88.7(13.4) 

IPSL 

25.0(2 1.6) 

66.4 (5.4) 

0.72 (0.80) 

47.1 (0.3) 

2 65.1 (0.0) 

2 66.4 (0.5) 

JMA 

27.0 (2 0.4) 

62.3 (4.9) 

0.31 (0.70) 

48.7 (7.2) 

2 122.8(2 8.4) 

2 122.4 (2 8.5) 

RACMO 

20.2 (2 3.5) 

68.2(11.9) 

0.40 (2 0.20) 

34.3(8.1) 

2 33.4(2 6.2) 

2 27.6 (2 6.2) 


"Numbers in the parentheses are the changes in the pertmbed climate. Listed are sensi ble and latent heat fluxes (SH, LH , in W/m2J, precipita- 
tion (PR EC, mm/day) , total cloud water path (TGLWP , g/m 2 ), shortwave, and total cloud-radiative effect (SWCRE , CRE, W/m 2 ). The asterisk 
denotes that the model has not reached equilibrium state (the ECETH model). 
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Table 5. Same as Table 4 but for Si I 


ModellD 

SH 

LH 

PR EC 

TGLWP 

SWCRE 

CRE 

ACCESS 

1 1.9(2 1.8) 

84.1 (7.4) 

0.26 (0.50) 

65.0 (2 10.8) 

2 123.0 (29.1) 

2 1 1 3.9 (26.4) 

CAM4 

23.7 (0.4) 

59.3 (7.9) 

0.00 (0.00) 

77.2 (4.8) 

2 133.4(2 1.7) 

2 129.7 (2 1.4) 

CAM S 

15.1 (2 0.3) 

90.2 (9.1) 

0.00 (0.20) 

55.0(14.9) 

2 124.1 (2.3) 

2 116.4(2.8) 

CCC 

29. 7 (2 4.4) 

63.3 (22.8) 

0.70 (2 0.33) 

228.2 (2 76.8) 

2 107.2 (14.8) 

2 100.8(17.4) 

CLUBB 

4.2 (0.7) 

88.5 (8.2) 

0.00' (0.00) 

25. 3 (6.3) 

2 95.7 (2 14.7) 

2 78.5 (2 13.6) 

ECHAM6 

2 21.4(1.7) 

78.4 (5.6) 

1.33 (0.90) 

173.1 (3.0) 

2 1 50.8 (0.4) 

2 1 50.9 (0.7) 

ECMWF 

6.8 (2 0.6) 

87.2(12.3) 

0.80 (0.13) 

48.7(15.1) 

2 24.6 (2 7.2) 

2 17.3(2 6.3) 

EC ETH 

6.5 (5.3) 

73.1 (15.4) 

0.31 (0.39) 

144.4 (35.0) 

2 129.4(2 7.1) 

2 130.1 (2 3.5) 

GFDL A M3 

15.5(2 6.3) 

78.7(15.8) 

0.30 (0.50) 

40.0 (5.5) 

2 118.4(2 11.8) 

2 111.3(2 11.2) 

GISS 

10.8 (0.5) 

76.3 (5.4) 

0.43 (2 0.07) 

129.8(2 55.7) 

2 68.1 (25.8) 

2 66.2 (25.1) 

GMAO 

5.1 (2 0.1) 

84.9 (8.0) 

0.70 (0.11) 

10.0(10.7) 

2 1 1.8(2 9.8) 

2 10.2(2 8.2) 

HadGEM2 

7.4 (2 0.2) 

69.7 (5.3) 

0.00 (2 0.00) 

4.7 (2 3.3) 

2 27.3 (9.4) 

2 24.5 (8.7) 

IPSL 

23.5 (2 0.9) 

74.7 (7.6) 

0.76 (0.1 1) 

52.2 (-0.0) 

2 69.3 (0.4) 

2 63.5 (0.7) 

JMA 

26.9 (2 4.1) 

73. I (14.5) 

0.59 (0.50) 

80.3 (20.1) 

2 157.1 (2 13.9) 

2 151.5 (2 15.4) 

RACMO 

15.3(2 2.3) 

91.0 (8.9) 

0.24(2 0.18) 

100.9 (2 5.5) 

2 85.7 (4.5) 

2 74.7(4.8) 


as an illustration of why the SCMs simulated different 
clouds in the control climate. 

5. Cloud Feedbacks 

5.1. SCM Results at SU (Stratocumulus) 

[21] We first use the cumulus under stratocumulus 
regime SI 1 to esta blish a framework to interpret the 
cloud feedbacks in the 15 SCMs. Figure 7 shows the 
change of net CRE from CTL to P2S at SI 1. Increase of 
CRE in the figure means positive cloud feed backs; 
decrease of CRE means negative feedbacks. For simplic- 
ity, the change of CRE is referred to as cloud feed back. 
The 1 5 SCMs simulated negative and positive cloud feed- 
backs that span a rather wide range of about 40 W/m 2 . 
Blossey et al [2013] showed this range as about 10 W/m 2 
in LES models. Because of the simplified CGILS setup, 
we do not expect the feedbacks here to be the same as in 
the full GCMs, but they allow us to gain some insight 
into the physical processes that detemune them. 

[2s] In Figure 7, the character "X" above a model's 
name indicates that shallow convection is not triggered 
in both the CTL and P2S sim ulations of t his model. 
The character "O" above a model's name indicates that 
shallow convection is active in at least one of the simu- 
lations of CTL and P2S. PBL schemes are always trig- 


gered in all mod els. Models without these characters 
about t hei r names used unified schemes of turbulence 
and shallow convection (such as CLUBB and 
RACMO) or did not subnu t informa tion for convection 
(such as ECMWF). One can see that models without 
active shallow convection tend to simulate negative 
cloud feed backs , while models with active convection 
tend to simulate positive cloud feedbacks. 

[29] Without convection, as discussed in the previous 
section for the JMA model, the water vapor balance is 
achieved by a competition between the moistening 
effect of the "turb" term in eq uation (7) and d rying 
effect of the net large-scale condensation "c-e" term and 
large-scale forcing; clouds are ca used by the moisteni ng 
term from the PBL scheme. Therefore, the response of 
the PBL scheme to SST largely determines the change 
of cloud water , hence, the cloud feed backs. Even 
tho ugh cloud microphysical and precipitation processes 
can also infl uence cloud feed backs, as mentioned 
before, since precipi tation is typically proportional to 
cloud wa ter, cloud wa ter controls the net change of con- 
densates in the simulations. 

[30] The PBL moistening term at the alti tude of maxi- 
m um cloud liq uid water is larger in the warmer clima te 
in virt ually all models as shown in Figu re 8a. In the one 
exception of the CCC model, the simulated altitude of 


Table 6. Same as Table 4 but for S6 


M odel_ID 

SH 

LH 

PREC 

TGLWP 

SWCRE 

CRE 

ACCESS 

6.8 (2 0.4) 

1 1 1 .4 (10.9) 

1.02(0.16) 

19.8 (0.9) 

2 9.6 (2 0.4) 

2 9.0 (2 0.4) 

CAM4 

8.5 (0.0) 

105.3 ( 12.2) 

0.00 (0.00) 

247.9 (24.0) 

2 177.4(2 4.5) 

2 160.1 (2 5.6) 

CAMS 

6.5 (2 0.2) 

104.3(13.4) 

0.74 (0.16) 

24.3 (2 3.4) 

2 35.2 (8.2) 

2 34.2 (8.1) 

CCC 

9.0 (0.5) 

122.4 (7.3) 

1.59 (0.60) 

68.9 (2 34.2) 

2 35.4 (24.9) 

2 27.3 (18.7) 

CLUBB 

10.4 (2 0.1) 

1 19.5 (10.2) 

0.57 (2 0.10) 

31.8 (2 0.6) 

2 91.7 (1.4) 

2 73.7 (2 0.1) 

ECHAM6 

2 5.6 (2 0.7) 

102.5 (9.2) 

0.79 (0.00) 

183.2(8.2) 

2 181.6(2 0.1) 

2 146.6 (2 4.1) 

ECMWF 

7.9 (0.6) 

108.1 (8.5) 

0.86 (0.70) 

25.5(6.0) 

2 12.6 (2 2.5) 

2 7.1 (2 2.5) 

EC ETH 

2.8 (2 1 .2) 

104.8 (7.6) 

0.61 (0.10) 

1 30.0 (5.0) 

2 125.5 (2.8) 

0.0 (0.0) 

GFDL A M3 

8.8 (2 0.6) 

1 10.0 (9.3) 

0.84 (0.12) 

5.9 (I.I) 

2 12.7(2 13.3) 

2 1 1.7 (2 13.4) 

GISS 

1 1.4(2 0.9) 

125.6 (1 0.0) 

1.41 (0.1 1) 

18.8(2 4.2) 

2 41.9 (12.8) 

2 39.4(11.4) 

GMAO 

6.1 (2 1.9) 

1 16.5 (6.1) 

1.14 (0.1 1) 

59.0(1.3) 

2 37.4 (1.0) 

2 33.1 (0.8) 

HadGEM2 

6.0 (2 0.4) 

109.9 (9.9) 

0.98 (0.12) 

3.7 (0.7) 

2 22.0 (2 0.9) 

2 20.1 (2 I.I) 

IPSL 

10.2 (2 0.5) 

1 1 8.8 (10.7) 

1.34 (0.17) 

74.6 (2 1.6) 

2 59.0 (4.7) 

2 53.6 (4.0) 

JMA 

14.7 (2 0.1) 

108.2 (7.8) 

0.63 (0.70) 

179.8(25.3) 

2 107.0 (2 5.9) 

2 101.2(2 6.3) 

RACMO 

12.0 (2 0.5) 

I 08.2 (8.1) 

0.66 (0.60) 

63.2 (7.6) 

2 28.4 (2 1 .6) 

2 25.8 (2 1.7) 
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Eigure 7. (a) Change of cloud-radiative effect (CRE, 

W/m 2 ) in SCMs at location SI I corresponding to 2 K 
SST pert urbation. Character M X " above a model’s name 
indicates that the shallow convection scheme is not 
active; "O" indicates that the shallow convection 
scheme is active . Models without these characters either 
do not separately parameterize shallow convection and 
PBL turbulence, or do not submi t results with convec- 
tion information. 


maxim um cloud wa ter in P2S is m uch higher than in 
CTL, above the top of the boundary layer (not shown) , 
where the turbulent term is small . The increased mois- 
tening by the PBL schemes is generall y consistent with 
the increase of surface latent heat flu x (LHF) in P2S, as 
shown in Figure 8b. The increase of latent heat fl ux 
wi th SST is consistent wi th CGILS LES sim ulations in 
Blossey et al . [2013] (their Table 4) and in earlier LES 
studies under similar experimental setup [e.g., Xu et al, 
2010]. AJso, Liepert and Pr evidi [2012] showed that in 
virtually all 21 st century climate change simulations by 
CMIP3 models, surface latent heat fluxes are larger in a 
warmer climate over the oceans (their Table 2, column 
3). 

[3 1] Previous studies [e.g., Caldwell and Bretherton , 
2009] have shown negati ve cloud feedbacks in mixed 
layer models (MLM) and have attribu ted the mecha- 
nism to la rger surface la tent heat fl ux and weaker large- 
scale subsidence i n a warmer climate. These two condi- 
tions are also shown in the CGILS SCM models that 
do not trigger convection . Table 5 shows that cloud 
water path in the negative feedback models is increased 
in the warmer cljmate. The example in Figure 6a 

(dashed lines), for the .TMA .model also, illustrates, the 
larger moistening rate by turbulence and deeper cloud 
layer in the warmer chnrnte. The CGILS results are 

therefore consistent wi th the interpretation of the nega- 
tive feed backs in MLMs . Exceptions are noted in wruch 
the convecti ve scheme is not acti ve in a model, but the 
model has small posi tive cloud feed backs, such as in 
CAMS and ECHAM6 . These may be related with 
cloud-top entrainment, included explicitly and im plic- 


Taking the ensem ble of models as a whole , we can use 
Figure 9a to schema tically summarize the negative 
cloud feed backs in the SCMs without convection. In 
these models, accompanied by the weaker large-scale 
subsidence, the warmer climate has greater surface 
latent heat flux, larger turbulence moisture convergence 
in the cloud layer, and consequently an inclination to 
give the negative cloud feedbacks. Tills mechanism is 
not new, but we see t hat it can explain the SCM resul ts 
in CGILS without activated convection. 

[32] We now turn to models wi th active shallow con- 
vection . Figure 7 shows that these models tend to have 
positive cloud feed backs . As discussed in the previous 
section for CAM4 and GISS, shallow convection acts 
to dry the cloud layer. It is a moisture sink that has the 
same sign as the stratiform condensation sink in equa- 


tion (7). The enhanced moisterung from the PBL 
scheme in the warmer climate is approximately bal- 
anced by enhanced drying from the sum of the strati- 
form condensation and shallow convection. If the rate 
of d rying from t he shall ow convection is greater than 
the rate of moistening from the PBL scheme as SST 
increases, the stra tiform condensation can decrease in a 
warmer clinmte. This tends to reduce cloud water and 
clouds, thus causing positive cloud feedback . The 
enhanced rate of convective drying in the warmer cli- 
mate may be explained by the moisture flux in equation 
(4) immediately above the top of the boundary layer. 
The moist ure contrast is larger in the wamler clima te, 
since the subsiding free tropospheric air remains dry 
but the total water in convective plumes increases with 
SST. An example is shown in Figure 6c for the GISS 
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p w w 

Model 

Eigur e 8. (a) Change of moisture tendency in the layer 
of maxim um cloud water {g/kg/day) by the " Turb” 
term from the control clinlate to the perturbed clima 
te at SI 1. (b) Same as Figure 8a but for surface 
latent heat 

flux (W/m 2 . ) 
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Control Climate Warmer Climate 



Negative 

Feedback 


(b) 




Positive 

Feedback 


(c) 




Feedback 

Uncertain 


Figure 9. Schematics of cloud feed backs. Changes of clouds from the (left) control to (right) warmer clima tes. 
Blue arrows denote the term of turbulence parameterization in the moisture budget eq uation; red arrows denote 
shallow convection. The sizes of arrows schematically correspond to the magni tude of moisture tendency from the 
associated processes . (a) Negative cloud feedback, dominated by the increase of surface turbulence, the "NESTS " 
negati ve cloud feed back mechanism (see text) . (b) Positive cloud feed back , domina ted by the increase of shallow 
convection or cloud-top entrainment, the "SCOPE " po siti ve cloud feedback mechanism (see text), (c) Cloud feed- 
back from shallow cumulus of sufficient depth, with sign depending on the cloud depth and lateral mixing. 


model by using the dashed lines. In the warmer climate, 
there is increase of t urbulence moistening, but larger 
increase of convective drying, and therefore red uced 
cloud water. Active convection therefore causes larger 
ventilation of the cloud layer in a warmer climate, 
which tends to decrease clouds and cause positive cloud 
feedbacks. This increase of convective mixing of bound- 
ary layer air together with the change of cloud -top 
entrainment ca uses more dil ution of the cloudy layer 
and therefore positive feedback. We can therefore use 
Figure 9b to schematically summarize the positive cloud 
feedbacks in the models. The net cloud feedbacks can 
be considered as due to two opposing roles of surface- 
based PBL t urbulence and shallow convection aided by 
cloud-top entrainment, wi th the latter dominating in 
most of the models in which convection is active. Figure 
9b also applies to models wi th parameterizations of sig- 
nificant cloud-top entrainment. The PBL scheme can 
also be dominant over the shallow convection scheme 
in some models, such as in CAM4. In this model, as dis- 
cussed in the previou s section, the peak drying of shal- 
low convection occurs below the cloud layer instead of 
within the cloud layer. 


[33] Brient and Bony [2012] used the larger moisture 
contrast between the free t roposphere and boundary 
layer in the warmer clima te to explain the positive 
cloud feedbacks in the IPSL SCM and GCM, while 
Kawai [2012] used the increased surface flux to 
explain the negative cloud feed back in the JMA SCM 
and GCM. These are consistent with the present inter- 
pretation. Figure 7 shows tha t in CGILS when con - 
vection is active, the positive feed back dominates the 
negative feedback. In GCMs or in the real atmos- 
pheres, any changes in the frequency of convection 
and convective mass fluxes would also matter. We call 
the a bove two competing mechan isms in Figure 9 as 
the "NESTS-SCOPE " (Nega tive feedback from Sur- 
face Tu rbulence under weaker Subsidence-Shallow 
Convection PositivE feedback) mechanisms . Obvi- 
ously, given the wide range of physical parameteriza- 
tions in models, this in terpretation may not fit all 
mod els. For exam pie, Zhang and Br etherton [2008] 
showed that in CAM3 the interaction of an unin- 
tended deep convection with the cloud microphysical 
scheme caused a negative cloud feedback in tha t 
model. Nevertheless, the delineation of the two 
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Eigure 10. Same as Figure 7, but for (a) S6, (b) SI 2. 
The models are ordered in the same sequence as in Fig- 
ure 7. One model (EC ECH) did not reach quasi- 
equilibrium state and it is ind icated by "NIA " . 


due to the com plications described above, convection at 
S6 does not necessarily correspond to posi tive cloud 
feedbacks. In au simulations, surface latent heat flux is 
greater in the warmer climate (Table 6). We may there- 
fore use the same framework as for SI 1 to think that 
the larger surface latent heat flux alone is a factor for 
more clouds in a warmer climate, but the other factors 
from shallow convection such as lateral mixing favor 
more dilution of clouds and a positive cloud feedback. 
The two effects compensate each other differently in the 
models because of the different assumptions in the spe- 
cific parameterizations. 

[36] Figure 1 Ob shows SCM results at S12, where SST 
is colder and subsidence is stronger than at SI 1. The 
corresponding changes of surface turbulent fluxes and 
cloud water path are given in Ta ble 5. Clouds a re 
restricted to wi thin the boundary layer. The sim ulated 
cloud feedbacks also span a wide range. Three models 

simulated no clouds at this location (GFDL AM3, .EC- 
ETH, and CAMS) (due to the constancy of forcing). 

Most models simulated the same cloud feedback signs 

as at SI 1. Some simulated opposite signs, one of which 


S6. Cloud feedbacks are generally consistent with the 
change of cloud liquid water path (Table 6). Partially 


competing mechanisms is a useful framework to inter- 
pret the majority of models. 

5.2. SCM Results at S6 (Shallow Cumulus) and at S12 
(Coastal Stratus) 

[34] We now use the same framework as we used for 
SI 1 to interpret SCMs results at the other two loca- 
tions. Before proceeding, we need to supplement our 
schematics with another scenario in which the depth of 
convection is large and mixing of cloud y ai r with dry air 
can occur la terally . If the cloud-scale dynamical fields 
and the environmental relative humidity are the same, 
larger d rying from convection is expected in P2S than 
CTL because of the larger difference of the absolute 
h umidity of moist ure across cloud lateral bounda ries 
just like across cloud tops. This is schematically shown 
in Figure 9c. Other factors such as cloud -scale dynam- 
ics, cloud depth , and cloud microphysics can also 
change in a warmer climate, leading to more compli- 
cated behavior of clo ud feed backs for t hicker clouds. 
This scena rio also incl udes regime change of clouds 
from stratocumul us to shallow cum ul us as exhibited by 
some models (e.g., CCC at SI 1, not shown). 

[3 s] Figure 10a shows the SCM cloud feedbacks at 
the shallow convection location S6, with a range of 
about 30 W/m 2 (in LES, models , the range is less than 3 
W/m\ The models a re ordered in the same seq uence as 
in Figure 7. Almost all model s sim ulated convection at 
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is the GISS model. As 
above the GISS model ii 
shallow convection is not active at SI 2, in contrast to 
be active at SI 1. Consistent wi th our hypothesis, the 
cloud feed back changed from positive to nega tive. 
The concept ual framework i n Figures 9a and 9b can 
be gen - erally applied to describe the behavior of cloud 
feed- backs in the SCMs at S12. 


J^LTS 


5.3. LES Results 

[37] The CGILS LES result s have been sum 
marized in Blossey et al [2013]. To compare wi th 
SCM results, in Figures lla-llc, we show the LES 
cloud feedbacks at the three locations of S6, S 11 , and SI 
2, respectively . The LES results are more consistent 
with each other than SCMs. A t the shallow cumulus 
location S6 (Figure I 1 a), LES model s simulated a small 
po sitive cloud feed- back except for DA LES and WRF 
that had negligi ble feedbacks. A t the stratocum ulus 
location SI 1 (Figure lib), all models except for SAM 
simulated positive cloud feedbacks. At the coastal stratus 
loca tion S12 (Figure 11 c), all except for DALES 
simulated negative cloud feed back. There is therefore 
consensus, but not uniform agreement, among the LES 
models with regard to simulated cloud feedbacks. 

[3s] Blossey et al. [2013] a ttributed the negative 
feed- back at S12 to the deepening of the cloud layer in 
a rela- tively well-mixed boundary layer that is related 
to weaker large-scale subsidence in the wamler 
climate. As mentioned before , this is also the 
interpreta tion of MLM negative cloud feedback and in 
the SCMs of CGILS as shown in Figure 6a. In some 
SCMs, vertical resolutions are not sufficient, so the d 
eepening of clouds cannot be Simula ted. In these 
models, the weaker subsi- dence leads to less su 
bsidence d rying in the warmer cli- ma te. This i s 
accom panied by larger turbulent convergence of moistu 
re into the cloud layer from enhanced surface flux and 
more liquid water. Therefore , the SCM interpretations 
are still consistent with the 
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Figure 11. Same as Figure 7 but in LES models, (a) 
S6,(b) SI I, and (c) SI 2. 


LES results of deepening boundary layer. At SI 1, Blos- 
sey et al. [20 13] attributed the positive feed back in the 
LES models to cloud thinning in a warmer climate 
ca used by decoupling of the boundary layer with the 
stratocumulus layer. In SCMs, the decoupled mixing is 
calc ulated by ei ther shallow convection or cloud-top 
entrainment or both, which has been shown to ca use 
positive cloud feed backs as in Figure 6c. At S6, Blossey 
et al [2013] attrib uted the posi ti ve feedback to more 
preci pitation. 

[39] A companion paper by Br etherton et al. [2013] 
investigated the sensi tivity of LES res ults to large-scale 
conditions, incl uding separate changes in su rface forc- 
ing, large-scale subsidence, environmental rela tive 
h umidity , and C0 2 concentration. These are not studied 
here since in CGILS we onl y aim at the total derivate of 
cloud feedback to imposed SST forcing with implied 
change in large-scale subsidence. The potential im pact 
of the chan ge of C0 2 forci ng is left for future stud y. We 
point out that the consensus among the LES models in 
Figure 11 does not necessarily mean they simulated the 
correct cloud feed backs. Nevertheless , they give pla usi- 
ble answers for SCMs to target for. Event ually, they 
need to be validated by observations under more realis- 
tic experimental setups. 


6. Summary and Discussion 

[4o] The experimental setup of CGILS was used to 
simulate shallow cumul us, stratocumulus, and coastal 
strat us and to investigate the physical mechanisms of 
cloud feedbacks under idealized climate change i n single 
column models. In models where shallow convection is 
not activated or plays minor role in drying the cloud 
layer, cloud feed backs tend to be negative . In models 
when convection is active, cloud feed backs tend to be 
posi tive in the stratocumu 1 us and coastal stratus regime , 
but uncertain in the shallow cum ulus regime. A frame- 
work is described to interpret the SCM cloud feedbacks 
by using the two opposing effects of increased moisten- 


ing from PBL scheme under weaker large-scale subsi- 
dence and enhanced drying from shallow convection in a 
warmer climate, with the former ca using negative cloud 
feed backs and the convective scheme ca using positive 
cloud feedbacks. The convective scheme plays a more 
dominant role at times when it is active. These mecha- 
nisms are sunlJDarized as the NESTS nega tive feed back 
and SCOPE positive feed back mechanisms. LES models 
simulated overall consistent positive cloud feed backs in 
the shallow cum ul us and stratocum ulus regimes, but 
negative feed backs in the coastal stratus regime. 
The LES results tend to support the NESTS-SCOPE 
mechanisms. 

[41] The relevance of CGILS results to cloud feed- 
backs in GCMs and in real-world climate changes is not 
clear yet. In a preliminary comparison to cloud feed- 
backs in four GCMs at the three locations, SCMs results 
were uncorrelated to those simulated by the parent 
GCM, suggesting the complexity of translating the 
results from SCMs to the feedbacks sim ulated by 
GCMs. While CGILS is motivated by understanding the 
physical mechanisms of cloud feedback s in GCMs, there 
are several issues that limit the applicability of the SCM 
results . First, the idealized forcing is steady state. Diur- 
nal and synoptic variabilities are not considered. Second, 
the large-scale fields are not interactive wi th clouds. 
Third, the spatial variability of GCM cloud feed back 
may be la rge and so direct comparison a t the selected 
locations may be inappropriate. Furthermore, the pat- 
tern of atmospheric large-scale condition in the GCMs 
may shift locations in a warmer climate [ Webb and Lock, 
2012]. Future phases of CGILS will investigate how 
results from the sinlplified case stud y should be used or 
how the case study should be modified to better under- 
stand cloud feedbacks in more complex model s and in 
observations. The CGILS results highlight the desira bil- 
ity to treat physical parameteriza tions in General Circu- 
lation Models (GCMs) as an integrated system rather 
than individual components in order to red uce cloud 
feedback uncertainties. 
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